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I. INTRODUCTION 


In a gaseous squeeze-film bearing, load capacity is generated by oscillating 
one of the bearing surfaces at a high frequency. An asymptotic analysis for 
large squeeze number was formulated in [1] and the performance for various 
squeeze-film bearing geometry was given in Refs. 2 to 5. It was shown that 
in order for the squeeze action to be effective the squeeze number must be 
large and the excursion amplitude be comparable to the mean film thickness 
In the literature cited above, the supported mass was assumed to be stationary 
disregarding the oscillatory forces imposed by the squeeze motion In reality 
the supported mass may respond dynamically with the squeeze action. This was 
investigated analytically in [6] by a simplified approach with the response 
assumed to be synchronous The results indicate that the load capacity may 
be considerably affected by the synchronous response depending on the squeeze 
frequency and the mass of the float. 

The supported mass motion was also studied in [71- For moderately high 
squeeze number, numerical solution for the response was obtained using small 
perturbation analysis The supported mass respond was found to be synchronous 
with the squeeze motion; this was also verified by experimental observation. 
The dynamic response of a spherical squeeze-film bearing with a time scale 
much larger than that of the squeeze motion was investigated in [8] using 
the results of [9]. 

The purpose of this report is to investigate the synchronous response of 
squeeze-film bearings with very high squeeze number. The geometry of a 
double squeeze-film circular disk has been chosen for simplicity. Three 
methods of solution are presented and the resul's compared. 
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II. ANALYSIS 


Consider a double squee ze- f i lm thrust plate as shown in Fig. 1. When the 
supported mass is at its central position, the mean film thickness on 
either side is C. Let l be the thickness of the thrust plate, then the 
mean distance between the two squeeze surfaces is 2C + i. Denote the 
instantaneous position of the lower squeeze surface, the thrust plate 
and the upper squeeze surface by y^, y 2 , and yy respectively. We have, 
then , 


yi = -C e cos co t (1) 

y 2 = C 6 ( 2 ) 

yj * 2C + 1 + C e cos to t (3) 


where e is the dimensionless excursion ratio and co, the squeeze frequency. 
Normally, if one neglects the dynamic response of the plate due to squeeze- 
film actions, the dimensionless displacement 6 would be time- independent . 
However, the dynamic response is, in general, not negligible. Therefore, 5 
is a periodic function of time. The time-average of 6 should be related to 
the steady-state displacement as follows: 


1 


2n 


2rr 


r 


C 6 d (ou t) = C - CT| 0 


(4) 


0 

where Ho is the dimensionless steady-state displacement from the central 
position; Tjo is defined to be positive for a downward displacement. 


The instantaneous film thicknesses of A and B are, from (1), (2) and (3), 


* y 2 ' yi = C 6 + C e cos a) t 


(5) 


fig = ■ y 2 ■ & - 2 C + C e cos u) t - C 6 


( 6 ) 


Assume that the squeeze number a is very high sc that we can apply the 
results of an asymptotic approximation. Let and Tg represent the 
asymptotic approximation of (ph) in film A and B respectively Using the 
results of [1], we obtain 


T A = Pa C 



(6 4- e cos 


(6 + e cos 


r) d t 


t) d t 


1/2 


(7) 


T B * P a C 


2tt 


. 2tt 


1/2 


(2 + e cos t - 6) d t 


f (2+6 cos T - 6) d T 


( 8 ) 


The instantaneous pressures in A and B can be easily obtained by 


Pa 


= *A = 


C (5 + e cos t) 


(9) 


P B = 



h 

C (2 +■ e cos t - 5) 


( 10 ) 


Note that although T A and Tg are t ime- i ndependent > p^ 3nd p B are functions 
of time. Furthermore, for a given excursion ratio e, the pressures p A a"d 
p D are not known, because 6 is an unknown quantity. Therefore, one more 

D 

equation is needed for 6, which is the equation of motion of the supported 
mass , 

. 2 

“7 TT2 * F y + P A A - P B A 


using Eqs. (7) to (10). Here, M is the mass, A is the bearing area, and F 
is the time-independent external force (positive upward) 


Define 


B = 


hL- 

p7^ 

M C ce 2 

~Pa A 


*A = 


f r2n 


1/2 


(6 + e cos t) d t 


■2tt 


- 


J 

0 


6 d t 


J 


( 11 ) 


X B 


C /-2n 


V/2 


(2 + e cos t - fj) d t 


r 


2tt 


TT 


6 d T 


V 0 

then, the equation of motion becomes 


B r^‘ ? y + 


6 + e cos t 


2 + e cos t 


( 12 ) 


integrating (12) over a complete squeeze cycle we obtain a relationship 
between Fy and e and 6- 


lA 

7n 



d t 

6 + e cos t 


+ 


ia 

2tt 



2 + e 


d t 

COS T - 6 


0 


( 13 ) 


With F y given by (13), Eq. (12) can be considered as the differential 
equation for If we are inteiested in the steady-state dynamic response 

rather than the transient, the boundary conditions to be satisfied are the 
periodicity conditions: 


6 (t) » 6 (t + 2 rr) 


04) 


d 6 (t) = d 6 (t + 2rr) (15) 

d t d t 


It should be pointed out here that Eq„ (12) is non-linear in 6* 
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III. METHODS OF SOLUTION 
A Harmonic Solution 

Assume that 6 c .1 be expressed by 

6 =■ 1 - T| 0 - X cos t - v sin T (16) 

where T] 0 is the steady-state downward displacement from the central 
position; X and v are the in-phase and out-of-phase components of the 
response. The higher harmonics have been truncated in (16). Substituting 
(16) into Eqs. (11), (12) and (13), we obtain 


= - Tl 0 ) 2 + | [(e - X) 2 + v 2 ]J 

- fd + v 2 + I O + *> 2 + ^ 2 ]J 


r (l - Tj 0 ) 2 - u - X) 2 - V 2 V(1 + • < e + X) 2 -v 2 j 


The equation of motion (12) takes the form 


B (Xcos t + v sin t) = F v + ~r, 1 ■ \ , — ; rr ; — 

y (1 - Tj 0 ) + (e - X) cos t v sin t 


(1 + T^) + (e + X) cos t + v sin t 


Multiply both sides of (18) by sin t and integrate from 0 to 2rr. 


/ sin t d r 

(1 - T) 0 ) + (e - X) cos t - v sin t 


ttB v * Ij 


sin t d t 

(1 + T] 0 ) + (e + X) cos t + v sin 
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Expanding the 
we obtain 


intergrands in power series of cos t and 

*-“Mo i - Uo 


sm t 


tt B 


v = 



v 3 v j 

' e - X 

2 3 

. 3 / V \ 

1 - Ho " 

1 ~ \ + 4 1 - Ho ( 

L 

i l - 

+ ‘I‘- T J + | 

lB -f 

V 2 V 

e + X’ 

2 i , 3 

l _ 2 M 1 n(n 

1 + ^io "L 

l + Hq ^ l + Tj 0 

U V 

1 4 J • 


Take the first power of v in the right hand side of (19) then 



Since the quantity in the bracket is, in general, not equal to zero, the 
only way to satisfy (19a) is V = 0. Furthermore, because the right hand 
side of (19) is a converging power series in V , we conclude that 


V = 0 


( 20 ) 


Using (20) and multiplying both sides of (18) by cos t and integrate from 
0 to 2rf, we obtain 


rr B X = I A 


I 


B 




- 2tt 


L±Jh 

e + X 


1 

YU - Ho > 2 - (e 



1 

Yd + Ho) 2 - 


(e + X) 2 




( 21 ) 


which is the equation to be solved for X Before we solve Eq. (21) 
which is obviously non-linear let us first obtain a linearized 
solution . 

For the case that 
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f-4 < < L 

1 - TV, 


and 


■ §~± - A - / < 1 
1 + Ho * 

we can linearize Eq. (21) and obtain 

2JW 

1-Tln 2 


02 ) 


X = 


(23) 


1 -Tlo 


T 


It is seen that when B = ~ — :ro , \ becomes infinitely large; the response of 

i- 

the mass is unbounded according to the linear approximation of a truncated 

harmonic analysis. Suppose that T ] 0 is positive, i.e., the mass has a steady- 

2 

state displacement downward. If B < ^ — , X has the same sign as e and 


1“T1 0 


the response is in-phase with the squeeze motion of the lower bearing surface. 

If B > ~ — —£ , the response is 180° out-of-phase with the squeeze motion and 

L ~ Mo 

the squeeze action is enhanced there. 

If condition (22) is not valid, we have to solve the nonlinear Equation (21). 

An iterative method may be used taking the linearized results as the initial 
guess for the iterative process. The numerical computation has been programmed 
on a computer. 

B. Exact Numerical Solution 


Rewrite Eq. (12) in the form. 


= G (6) 


d t 


where G ( 6 ) = B ' 1 


/ 


| F y + 


IA 


'B 


2 4- e cos t - 5 


(24) 

(25) 


6 + e cos t 
Fy, I A and I g are given by (11) and (13). 

An iteration method will be used to solve Eq . (24). Let be the i-th 
approximation for the solution of (24). Let [ be a correction function to 
6 i, so that the sum, 

6 i+l = §i + u i (l 


Si + Ui 
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will be an improved approximation Thus 


6i d 2 uj 

d t 2 d 


- G (6j + Uj) 


(26a) 


Taking the first two terms of Taylor's Series expansion of the right hand 
side, we have 


d 2 £,- + d 2 Ui 

d t 2 d t2 


= G (6t) + u t G' (6i) 


(27) 


or 


- UjL G' (6i) = G (6,) - *Lii 


d 


d t' 


(? 7 a) 


Since 5j is a periodic function and so is 6 j + j : Eq . (26) then indicates that 
uj must also be periodic Thus 


Ui (t) = Uj ( t + 2tt) 
d ( T ) - d m ( T -K 2 tt> 

d * d T 


( 2 *> 


Note that G (6i) is a function of 6i as well as a functional of 5i because 
6i is implicitly involved in the definite integrals I A and Ig For the 
iterative process, an approximate expression for G'(6i) can be obtained bv 
treating I A , Ig and Fy as independent of 5 X Thus, 




(2 + e cos t + 6i)‘ 


( 29 ) 


Since the approximate expression (29) is used to determine Uj for the nex<- 
approximation in 6i, it would affect the speed of convergence of the iterative, 
process; however, it would not affect the accuracy of the resulting final 
solution . 


Let there be Q points between 0 and 2rr, then the interval is 


A = 


2tt 


(30) 


Q - 1 
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In finite-difference form, 



The superscripts q, q + 1 and q - 1 indicate stations Substitute into 
Eq- (27), 


_ q q - 1 
2 u. M + u^ 


= A 2 G (6i q ) - [6i q + 1 - 2 6i q + 6i q ' L ]; q = 1,2,..Q-1 


Or, in matrix form 
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V = 


A^ G (6j ) - [6^ - 2 6 £ + 6^ J 

A 2 G (Si 2 ) - [6i 3 - 2 Si 2 + Si 1 - 

f 

I 

A 2 G (6i q “ 2 ) - [S^" 1 - 2 6i q " 2 + 6i q ' 3 ] 

A 2 G (S^" 1 ) - [6 ± q - 2 S i q " 1 + 6i q " 2 ] 


(33) 


We have used the peroidicity conditions 








(34) 


(35) 


To start the iteration, the result of the harmonic solution, (16) can be 
used as an initial guess. The u^ ’ s can be easily calculated from (31) by 
inverting the matrix M. Repeat the process as indicated in (26) until the 
elements of Ui's are less than a specified value. 


Typical results are obtained for a squeeze frequency of 18000 rad/sec, 
e=0.3, Tj 0 = 0.2. C = 0.001 in, A = 1 in^ , p a = 14.7 psi and M = 0 1 lb. 
The 2tt interval is equally divided into 25 parts (Q = 26). The results of 
the numerical solution are tabulated as follows: 
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t 

0.0000 
0 2513 
0.5027 
0.7540 
1 .0053 
1.2566 
1.5079 
1.75S2 
2 .0106 
2 .2619 
2 5132 

2 7646 

3 0159 
3.2672 
3.5185 
3.7699 
4.0212 

4 2725 

4 5238 
4. 7752 

5 0265 
5 2778 
5 5291 
5.7805 
b .03 18 
6.2831 


6 (t) 

0.88907 
0 8869^ 
0.88058 
0.87012 
0.85578 
0.83793 
0.81711 
0.79412 
0.77010 
0.74662 
0.72571 
0.70975 
0 70104 
0.70104 
0.70975 
0. 72571 
0.74662 
0.77010 
0-79412 
0 81711 
0 83793 
0.85578 
0.87012 
0 88058 
0.88694 
0 88907 


Thus, from Eq. (16) 

6(t) = 1 - T] - X cos t 
1 o 

= 0.8 + 0.0884 cos t 

On the other hand, the non-linear harmonic solution yields X = -0.0884. 

The linearized harmonic solution yields an X of -0 035- Ic can be easily 
seen that only the non-linear harmonic solution agrees well with the 
tabulated numerical solution whereas the linear harmonic solution does not 

LOAD CAPACITY 


We can identify (-Fy) as the dimensionless load capacity W of the double 
squeeze-film thrust plate. Several expressions for W can be obtained 
depending on how the synchronous response is calculated; they are 
enumerated as follows; 

1) Load capacity without taking the synchronous response into account 
can be obtained from Eq. (17) by setting X =0. 
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W =-F = 

O r yo “ 


- V 2 + 7 e 2 

a - v 2 - 7 





+ 7 


(1 + v 2 - * 2 


(36) 


2) Load capacity using the harmonic solution for the synchronous response 
is directly obtainable from Eq . (17). 


W = F y = 


(I - Tip) 2 + y (e - X) 2 
(1 - Tl 0 ) 2 - (e - \) 2 


((1 + T| n ) 2 4- 7 ( e + X ) 2 

(i + v 2 ' <e + X ) 2 


(37) 


The load capacity is designated by and for respectively the linear and 
non-linear results of 

3) Load capacity Wg can be calculated from Eq. (13) using the exact 
numerical solution for 6 - 


W E = r 
2tt 


2 rr 

J 77 


d T 


e cos t 


IB 

2tt 


. 2 rr 


d t 


O' 


2 + e cos t 


(38) 


In order to visualize the effect of synchronous response of the supported 
mass to the load capacity of a squeeze-film bearing, it is illustrative to 
normalize the load capacity with respect to W Q . In Figure 2. Wg/W Q , Wjj/W 0 

and Wl/W q are plotted against B for e = 0.3 and 7 | 0 = 0 2 It is seen that 

when B is large, the three normalized load capacities approach to unity 
asymptotically. This indicates that for large B, the amplitude of the 
synchronous response is small so that its effect to load capacity is negli- 
gibly small. It is interesting to note that in the region where B > 6 , 
the normalized load capacities are positive and greater than unity; this 
is because the respcnse is 180° out. of phase with the motion of the closet 
squeeze surface (the lower surface if TIq is positive). This out-of-phase 
response enhances the squeeze action of the lower film A, while decreases 
the squeeze action of the upper film B Because the lower film is dominating 
for positive Tl 0 , the net effect is an increase in load capacity However, 
if B is small (B < 2, say), the response is in phase with the lower surface 

which results in a decrease in load capacity. As seen in Fig 2, the 

normalized load capacity becomes negative which is clearly undesirable 
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The response of the supported mass is seen to be similar to the response 
of the forced motion of an undamped spiing-mass system. When the frequency 
is lower than the critical frequency the response is in-phase with the 
forcing function. But, when the frequency is higher than the critical 
frequency (super critical), the response is 180° out-of-phase with the forcing 
function It is in this supercritical region we should design and operate a 
squeeze-film bearing. 

In Fig. 2, it is clear that is very close to Wg, especially in the super 
critical region. In other words, the non-linear harmonic solution agrees 
well with the exact numerical solution. This feature will be utilized in 
future squeeze-film bearing analysis (Synchronous Response of Conical and 
Spherical Squeeze-film Bearings). 

In Figures 3 and 4, W£/W Q is plotted for e = 0.1 and 0 3 respectively with T| 0 
as a parameter Note that in the supercritical region the normalized load 
capacity increases with T| Q . This indicates that the synchronous response 
increases not only the load capacity but also the stiffness. 
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IV. DISCUSSIONS 

In the exact numerical solution of the dynamic response of a double-film 
squeeze bearing, periodic boundary conditions for the function and its first 
derivative have been assumed at t = 0. In so doing, we have included the 
synchronous response and its higher harmonic responses, whereas the sub- 
harmonic responses are automatically discarded 

In order to investigate the importance of the sub-harmonic responses it is 
only necessary to impose the periodic boundary conditions at 0 and n2rr 
Clearly, for n = 2, for example, the solution would include half-harmonic, 
synchronous and higher harmonic responses. 

An actual plot is shown in Fig. 5 for e=0 3 and Tj Q - 2; the solid curve is 
for n = 1 and the dotted one for n = 2. Note that although the periodic 
boundary conditions are imposed at 0 and 4it for the n = 2 case, the actual 

numerical solution shows that it repeats itself from 2rr on. Also, the fact 

that the two curves being so close together indicates that the response is 
mainly synchronous. The computed load capacity, Wg/WQ, is 2 831 for n = 1 
and 2.906 for n = 2. Therefore, it is not necessary to take the sub -harmonic 
responses into consideration. 

Another point of significance is that the pressure forces in a squeeze film 
were ODtained by an asymptotic approximation In so doing, the damping part 
of the squeeze-film force becomes zero. The damping can be realized by 
considering the edge correction of the asymptotic approximation as will be 
seen in the following. Starting with the Reynolds' equation of a squeeze 
film thrust plate. 

H 3 J&_. (pr Ml = a .Zim (39) 

R BR \ dR 9 t 

where P, R and H are normalized quantities 


Let \|f = PH 
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rtn asymptotic analysis shows that we can write 

\|r (R, t) = ^00 (R) + * e (R, t) (41) 

Here, ^ (R) is the asymptotic approximation for large o, and ^ (R, t) 

is the edge correction which is important only in a region near the edge; 
the extent of which is of the order of 1/ V7 

Ihe governing equation for is 

2 2 

2? ^ = typ (42) 

2 d § 2 a T 

with boundary conditions 

♦e (? = °> t) = H 

t) = ^ ( R = 1) (43) 

and the periodicity condition 

♦e (') = (t + 2rr) (44) 

where § - (1-R) (45) 

From the structure of (44), it is not difficult to see that if H involves a 
cos t variation in time, would have both sin t and cos t terms. 

Therefore, when we obtain the instantaneous film pressure from 

p (R, T ) = - Y <v (R) + ¥ e (46) 

H 


it is clear that \|i e (R, r) would contribute to a damping. Thus, the 
damping occurs only in the edge region and is of the order of 1/ V7 



The stability of this double squeeze-film thrust plate [as been considered 
using the same approach ! .s ^10],i.e. the small parameter .,1 bility technique. 
Starting from Eq . (12), it is found that (See Appendix I) for e = 0.3 
and T] = 0.2, the bearing would be stable if B 2 12. As shown in Fig. 2, 

B = 12 is already in the supercritical region. Actually, this stability 
result is conservative because damping has been neglected. 



V. CONCLUSIONS 


The dynamic response of a double squeeze-film thrust plate has been 
analyzed theoretically. Based on the results obtained the following 
conclusions can be drawn: 

1. The dynamic response is mainly synchronous; sub-'narmonic 
responses are negligibly small in magnitude. 

2. Damping of a squeeze-film occurs only in the edge region 

and is of the order of 1 /~\fo For large squeeze number, 

the damping can be neglected. 

3. When the squeeze frequency is high (large B) , the response 
is 180° out-of-phase with the squeeze action. This 
enhances the squeeze action, and thus increases the load 
capacity and stiffness. The converse is true for low 
squeeze frequency ('•mall B) . 

4. The harmonic solution for the dynamic response without 
linearization agrees well with the exact numerical 
solution The linearized harmonic solution has only 
qualitative agreement with the exact numerical solution. 
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APPENDIX I STABILITY CONSIDERATION 

The stability of the double squeeze-film thrust plate can he formulated using 
Eq. (12) and setting 

6 = 6 p + 6' (1.1) 

where 5 q = 1 - T| and 6 1 is a small perturbation from the equilibrium 
position 6 • 

Substituting (I i) into Eq (12) we have 


~I A 

’ 1 

e :os t| . Ig 

1 e cos T f| 

Jo 

Jo “ 2 

C) 2 - 6 o 

I 2 - «o < 2 - 6 o) 2 jJ 


( 12 ) 


Note that in the above equation, we have neglected non-homogeneous terms 

such as — (1 - cos t) and non-linear terms involving 5' , etc. This 

"o “° 

is justified for a small parameter stability analysis. 


Rearranging Eq. (1.2) Eq . (1.2) we can write 


where 


^ — £l. + (a^ + a 2 cos t) 6' = 0 
d t *~ 



(1.3) 


(I 4) 


Clearly, the sign of a 2 is irrelevant because a 2 cos T is only a 
simusoidally fluctuating term. Neglecting the squares of 6' I A and Ig 


can be written as 
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Ia = 6, 


3 2 

+ T e z 


2 3 2 

*B = < 2 " V if / 


Thus , 


3 1 = 


H 

e 

2 3 
+ 1+ 7 

, I 2 ] 



1 7] 

o 




r « _ 

('* 

3 / € 

l 2 ] + i fl + 3 

e 


o 

l 

*—4 j 

2 I 1 -’cl 

1 J 1 + \ { 2 

1 + nJ 

iJ 


(1.5) 


( 1 . 6 ) 


Note that we have used the asymptotic approximation for large q. In so 
doing, the damping in the squeeze-film is neglected. The stability 
result is therefore conservative. 

Results of [10] and [11] indicated that for marginal stability 

a l = 0.23 - 0.5 a 2 (I-7> 

If k = 0.3 and TJ o =0 2, Eq. (1.7) yields B = 12. Therefore, for stable 
operation B should be greater than or at least equal to 12. 
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NOMENCLATURE 


A 

B 

C 

e 

F y 
y 

G 

h A’ h 3 

V *B 

M 

P 

Pa 

P 

r 

r o 

R 

T A> T b 
v 

W 

u 

"o 

w, w w 

h L’ W N* W E 


yi» y 2 > y 3 


6 

A 

e 


X, v 

P 


area 

MCco 2 /p a A 

mean film thickness when bearing is unloaded 

amplitude of excursion 

external force 

F/d A 
> a 

defined in (25) 

film thicknesses given in (5) and (6) 
defined in (11) 

mass of float: matrix defined in (32) 
pressure 

ambient pressure 

P/P a 

radial coordinate 
radius of plate 

r/r o 

asymptotic approximations of ph 
defined in (33) 
dimensionless load capacity 

" " " without supported mass motion 

11 " capacities with supported mass motion 

calculated by linear and non-linear harmonies 
solution and by exact numerical solution 
displacements defined in Fig. 1 
and Eqs. (1), (2) and (3). 

dimensionless displacement of supported mass 

defined in (30) 

e/C 

dimensionless steady-state displacement of supported 
mass from central position 
amplitudes of mass motion 
viscosity 


squeeze number 


u)t 

PH 

squeeze frequency 
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